MAPPING PROPERTIES OF THE HEAT OPERATOR ON EDGE 

MANIFOLDS 



ERIC BAHUAUD, EMILY B. DRYDEN, AND BORIS VERTMAN 

Abstract. We consider the heat operator on spaces with complete and incomplete 
edge metrics. We establish its mapping properties, recovering and extending the 
classical results from smooth manifolds and conical spaces. The estimates yield short- 
time existence of solutions to certain semilinear parabolic equations. Our discussion 
reviews and generalizes the earlier work by Jeffres and Loya. 
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1. Introduction and statement of the main results 

Edge metrics of complete and incomplete type provide interesting examples of ge- 
ometries that include asymptotically hyperbolic, asymptotically cylindrical and conic 
spaces. Incomplete and complete edge metrics behave differently from a geometric 
perspective. However, these metrics lend themselves to a parallel approach when one 
wants to do constructions of an analytic flavour. For example, to analyze asymptotics 
of solutions to the heat equation near spatial infinity in the complete case, one must 
introduce a geometric compactification. The resulting complete edge metrics, like their 
incomplete cousins, give rise to differential operators that are degenerate or singular in 
a sufficiently controlled manner that generalization of classical results is possible. In 
particular, we consider the Laplacian on manifolds with either type of edge metric and 
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establish regularity properties of solutions to the heat equation. We study the cases 
in parallel, reviewing and extending previous work by Jeffres and Loya for conic and 
b-metrics [JL03, JL04]. _ 

We begin by defining these metrics precisely. Let M be an m-dimensional compact 
manifold with boundary dM, where dM is the total space of a fibration <fi : dM — > B, 
and the fibre F and base B are closed manifolds. We consider a defining function 
x : C°°(M) R+ U {0} of the boundary dM with ar^O) = dM and dx ^ on dM. 
Using the integral curves of grad(x) we identify a collar neighbourhood (JcMof the 
boundary dM with [0, 1) x dM, where dM is identified with {0} x dM. 

Definition 1.1. A Riemannian manifold (M\dM, g) := (M,g) has a simple non- 
iterated 

(i) incomplete edge at B if, with respect to the identification U\dM = (0, 1) x dM, 
the Riemannian metric g = go + h and go attains the form 

g \ U\dM = dx 2 + x 2 g F + (j)*g B , 

(ii) complete edge at B if, with respect to the identification U\dM = (0, 1) x dM, 
the Riemannian metric g = go + h and go attains the form 

dx 2 6*q B _ 
go U\dM = — + ^ r + g F , 
x 1 x z 

where g B is a Riemannian metric on the closed manifold B, g F is a symmetric 2- 
tensor on the fibration dM restricting to a Riemannian metric on each fibre F, and 
\h\ go = 0(x) as x —7- 0. 

We will refer to either of these metrics as an edge metric. As the names imply, in- 
complete or complete edge metrics are incomplete or complete as Riemannian metrics. 
In the complete case dM is at infinite distance from any point in the interior of M; 
note that the two types of edge metrics are conformal in the interior of M. Exam- 
ples of complete edge metrics include asymptotically hyperbolic (conformally compact 
or 0-metrics) and asymptotically cylindrical (b-metrics) metrics, with dimF = and 
dim B = respectively. The product metric on H fe+1 x S-^ provides an example of a 
complete edge metric with trivial fibration structure. Examples of incomplete edge met- 
rics include conic metrics and conformal compactifications of asymptotically hyperbolic 
metrics. 

We consider a slightly restricted class of edge metrics, requiring that (ft : (dM, g F + 
4>*g B ) — > (B, g B ) be a Riemannian submersion. Recall that if p G dM, then T p dM 
splits into vertical and horizontal subspaces as dM © T^dM, where T^dM is the 
tangent space to the fibre of </> through p and T^dM is the orthogonal complement of 
this subspace. The new condition on go implies that the restriction of the tensor g F to 
T^dM vanishes. Moreover, in the case of incomplete edge metrics we need to assume 
that the Laplacians associated to g F at each b e B are isospectral. We summarize these 
additional conditions in the definition below. 

Definition 1.2. Let (M, g) be a Riemannian manifold with an edge metric. This metric 
9 — 9o + h is said to be feasible if 
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(i) (f) : (dM,g F + <p*g B ) — > (B,g B ) is a Riemannian submersion; 

(ii) if the edge metric is incomplete, the Laplacians associated to g F at each b E B 
are isospectral. 

Our perspective for establishing regularity for solutions to the heat equation is as 
follows: the mapping properties of the corresponding heat operator are encoded in the 
asymptotic behaviour of the heat kernel on an appropriate blowup of the heat space, 
as established by Mazzeo and the third author [MV10] for incomplete and by Albin 
[Alb07] for complete edge metrics. The regularity is discussed in terms of spaces with 
bounded edge derivatives and appropriate Holder spaces, which take into account the 
underlying singular or asymptotic geometry. To make this precise we introduce the 
notion of edge vector fields V e . 

Definition 1.3. Let M be a compact manifold with boundary dM being the total 
space of a fibration <fi : dM — > B with fibre F. Define the vector space V e to be the 
space of vector fields smooth in the interior of M and tangent at the boundary dM to 
the fibres of the fibration. 

The space V e is closed under the ordinary Lie bracket of vector fields, hence defines a 
Lie algebra. In local coordinates V e can be described as follows. Let y — (yi, y b ), b — 
dimi? be the local coordinates on B lifted to dM and then extended inwards. Let 
z = (zi, Zf), f = dimF restrict to local coordinates on F along each fibre of dM. 
Then (x, y, z) are the local coordinates on M near the boundary and the edge vector 
fields V e are locally generated by 

d d d d d 



dx dy± dyb dz\ dzj 

Definition 1.4. Let {M,g) be a Riemannian manifold with an edge metric. 

(i) Let C° (M, g) denote the space of functions that are continuous up to the bound- 
ary of M. The space C*(M,g) with k E N consists of all u E C°(M,g) 
such that for any choice of edge vector fields Vj E V e , the edge derivatives 
V\ ■ ■ - V n u E C®(M,g) are again continuous up the boundary for any n < k. 
We put 

||w||fc := sup \u\ + sup \V\u\ + • • • + sup \V\ • • • Vuu\ < oo. 

MM M 

(ii) The Holder space C"(M, g),a E (0, 1), consists of continuous functions u such 
that 

ii ii ii, Kg)-^MI / ^ 

\\u L := sup \u\ + sup r < oo, 

m m d(p,q) a 

where d(p, q) represents the distance between p,q E M with respect to the 
Riemannian metric g. Note that in the local neighbourhood of the edge these 
distances are uniformly equivalent to 

d((x, y, z), (x, y, z)) 2 ~ \x — x\ 2 + \y — y\ 2 + (x + x) 2 \z — z\ 2 , 
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if g is incomplete, and 

7/ / \ / — ~\\2 ^ \ x ~ x \ \y ~ y\ , i ~i2 

d({x,y,z),(x,y, z)) « ? — + * - z , 

(x + x) z 

if g is complete. 

In the case of incomplete edge metrics it also makes sense to restrict the Banach 
space of continuous functions to those which are fibrewise constant at x = 0. This is 
precisely the space of continuous functions ^ e °(M, g) with respect to the topology on M 
induced by the Riemannian metric g. The corresponding space of continuous fc-times 
edge-differentiable functions shall be denoted by ^(M,g), with 

tf k (M, g) := {u G C k (M, g)\V 1 --- V jU G ^ e °(M, g) for any V, G V e }, 

which is a Banach subspace of C k (M,g). For the Holder space with fractional differen- 
tiability we set ^ e Q (M, g) := C"(M,g). Note that similar considerations are possible in 
the setup of complete edge metrics, but in that case these do not lead to a refinement 
of the regularity statement. 

Our estimates with respect to these spaces are as follows, where a precise definition 
of the heat operator is given in §2. 

Theorem 1.5. Let (M,g) be an m- dimensional Riemannian manifold with a feasible 
incomplete non-iterated edge metric. Let e~ tA denote the heat operator corresponding to 
the Friedrichs extension A of the associated Laplacian on (M,g). Then e~ tA is bounded 
for any fixed t G 1R + with estimates of the form 

(1.0.1) e"' A : C k e (M,g) C k e +n (M,g), \\e~ tA f\\ c ^ < Cr n / 2 \\f\\ cl , 

(1.0.2) e"* A : C?(M,g) C e 2 (M, g), \\e- tA f\\ cl < Cr^\\f\\ G? . 

Precisely the same estimates hold if C*(M, g) is replaced by ^{M^g). 

For complete edges the result is similar, although the heat operator acts between 
weighted edge spaces in this case. A function u is in a weighted edge space, denoted 
u G x w C k (M, g), if and only if u = x w v, with v G C k (M,g). 

Theorem 1.6. Let (M,g) be an m- dimensional Riemannian manifold with a feasible 
complete non-iterated edge metric. Let e~ tA denote the heat operator of the unique self- 
adjoint extension A of the associated Laplacian on (M,g). Then for any w G K, e~ tA 
is bounded for any fixed t G M + with estimates of the form 

(1.0.3) e~ tA : x w C k (M,g) -> x w C k+n (M, g), ||e-* A /||^ + n < CT^H/H^c*, 

(1.0.4) e~ tA : x w C?(M,g) -> x w C 2 e (M, g), \\e~ tA f\\ xWCl < Cr 1+a / 2 \\f\\ x ^. 

An immediate observation from these estimates is that after taking into account the 
geometry of the underlying space, the mapping properties of the heat operator resemble 
the well-known behaviour on compact manifolds. As a particular consequence of our 
results, we establish short-time existence for solutions to certain semilinear parabolic 
equations on manifolds with edge metrics. Applications of such equations, including 
the reaction-diffusion equation, may be found in [CH98]. 
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Corollary 1.7. Let (M,g) be an m- dimensional Riemannian manifold with a feasible 
edge metric. Let A denote the Friedrichs extension of the Laplacian on (M, g) for 
an incomplete edge metric, and the unique self-adjoint extension of the Laplacian for 
a complete edge metric. Let the pair (X,Y) denote either (C|(M, g),C^(M, g)) for 
a G (0, 1) or 

(i) ( < ^f +1 (M, g),^(M, g)), k G N, if g is an incomplete edge metric; 

(ii) (Cg +1 (M, g),C*(M, g)),k eN, if g is a complete edge metric. 

For u G X , let Q : X —>Y denote a locally Lipschitz inhomogeneous term that depends 
on at most one edge derivative of u. Then the initial value problem 

Ou 

— + Au + Q(u,Vu) = 0, u (0) = feY, 
at 

admits a solution u G C([0,T]; X) on some time interval [0,T] for T > 0. 

Our discussion of the heat operator in the edge setup reviews and generalizes the 
work of Jeffres and Loya in [JL03] and [JL04], where the case dim I? = was con- 
sidered. Their work was based on the heat kernel analysis by Mooers [Moo96] in 
the incomplete (conical) case and by Melrose [Mel93] in the complete (b-cylindrical) 
setup. The analysis of the heat operator in the presence of incomplete singularities 
was initiated by Cheeger [Che83], with major contributions by Briining and Seeley 
[BS87], [BS91], Lesch [Les97], Melrose [Mel93] and Mazzeo [Maz91], to select a 
few. Related questions on regularity properties of solutions to parabolic equations in 
the singular setup have been studied in [CSS02], [Li07] and [LoyOI]. Most recently, 
a study of the inhomogeneous Cauchy problem on manifolds with incomplete conical 
metrics has been presented by Behrndt [BehII]. 

The paper is organized as follows. In §2 we review the asymptotic properties of the 
heat kernel as a polyhomogeneous distribution on the appropriate blowup of the heat 
space. We discuss the incomplete and complete cases separately since the asymptotic 
properties are different. In §3 we apply the asymptotics of the heat kernel to derive 
the mapping properties of the heat operator and hence the regularity of solutions to 
the heat equation, carefully estimating the corresponding integral in various regions of 
the heat-space blowup. Finally, in §4 we explain how these mapping properties yield 
short-time existence for solutions to certain semilinear parabolic equations. 
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support by the German Research Foundation DFG as well as by the Hausdorff Center 
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2. ASYMPTOTICS OF THE HEAT KERNEL ON EDGE MANIFOLDS 

Recall that we let (M m , g) be a Riemannian manifold with an edge at B b and a feasible 
edge metric g, where M is the interior of a compact manifold M with boundary. This 
boundary dM is the total space of a fibration : dM — > B, with fibre F* . The local 
coordinates on M near dM are given by (x,y,z), where (y) are local coordinates on 
B lifted to dM and then extended inwards, and (z) restrict to local coordinates on F 
along each fibre of dM. Note that m = 1 + b + f. 

Denote by A the Friedrichs extension of the Laplacian on (M, g) in the case of an 
incomplete edge, and the unique self-adjoint extension of the Laplacian in the case of 
a complete edge. The heat operator of A, denoted e~ tA , solves the homogeneous heat 
problem 

f (d t + A g )u(x,y,z,t) = 

1 u(x,y,z,0) =f(x,y,z), 

with u = e~ tA f. It is an integral operator 

(2.0.1) e~ tA f(p)= [ H(t,p,p)f(p)dvo\ s (p), 

J M 

2 

with the heat kernel H being a distribution on the heat space Ml := M + x M ; local 
coordinates in the singular neighbourhood are given by (t, (x, y, z), (x, y, z)), where 
(x,y, z) and (x, y, z) are coordinates on the two copies of M near the edge. 

The heat kernel H has non-uniform behaviour at the diagonal D and at the subman- 
ifold 

f {(£, (0, y, z), (0, y, z)) G R + x (dM) 2 : t = 0, y = y}, if M is incomplete edge, 
\ Ht, (0, y, z), (0, y, z)) G R + x (dM) 2 :y = y}, if M is complete edge. 

The asymptotic behaviour of H near these submanifolds of M\ depends on the angle of 
approach to these submanifolds, and is crucial for estimation of the integral (2.0.1). This 
dependence on the angle in the asymptotics of the heat kernel is conveniently treated 
by introducing polar coordinates around A and D. Geometrically this corresponds to 
appropriate blowups of the heat space M%, distinct for the incomplete and complete 
cases, such that the corresponding heat kernel lifts to a polyhomogeneous distribution 
on that blowup space, in the sense of the following definition. 

Definition 2.1. Let 22J be a manifold with corners, with all boundary faces embed- 
ded, and {(Hi, Pi)}f =x an enumeration of its boundaries and the corresponding defining 
functions. For any multi-index b = (b\, . . . , b^) G we write p b = p b ^ . . . p h ^ . Denote 
by H(2U) the space of smooth vector fields on 2U which lie tangent to all boundary 
faces. A distribution w on 2U is said to be conormal if w G p b L 00 (2U) for some b G 
and V 1 ...V i w G p fe L°°(2IT) for all V d G V b (W) and for every £ > 0. An index set 
Ei = {(j,p)} C C x N satisfies the following hypotheses: 

(i) Re (7) accumulates only at +00; 

(ii) for each 7 there exists P 7 G N such that (j,p) G Ei for p < P 7 < 00; 

(iii) if (7,^) G E h then (7 + j,p') G Ei for all j G N and < p' < p. 
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An index family E = (Ei, . . . ,En) is an iV-tuple of index sets. Finally, we say that 
a conormal distribution w is polyhomogeneous on Q2J with index family E, denoted 
w G ^f g (2U), if w is conormal and if near each Hi we have 

W ~ a 7:PP/( 1 °S/ ? i) P ' as p, ; ->■ 0, 

with coefficients a 7 p conormal on Hi and polyhomogeneous with index Ej at any H^Hj. 

For more on polyhomogeneous distributions and other relevant background material, 
we refer the reader to the classical references [Mel93] and [Maz91], as well as the 
excellent introduction to the 6-calculus by Grieser [Gr99]. 

2.1. Heat kernel on spaces with incomplete edge metrics. To obtain the correct 
blowup of M% in the case of an incomplete edge metric, one first does a parabolic blowup 
of the submanifold A defined above. The resulting heat space [M|,A] is defined as 
the union of M%\A with the interior spherical normal bundle of A in M\. The blowup 
[M%, A] is endowed with the unique minimal differential structure with respect to which 
smooth functions in the interior of M\ and polar coordinates on M\ around A are 
smooth. This blowup introduces four new boundary hypersurfaces, which we denote 
by ff (the front face), rf (the right face), If (the left face) and tf (the temporal face). 

The actual heat-space blowup is obtained by a parabolic blowup of [M|, A] along 
the diagonal 

D := {t = 0, x = x, y = y, z = z} C M\, 

lifted to a submanifold of [M^,A]. The resulting blowup is defined as before by 
cutting out the submanifold and replacing it with its spherical normal bundle; a new 
boundary hypersurface, the temporal diagonal (td), appears. The blowup is a 
manifold with boundaries and corners as depicted in Figure 1. 




Figure 1. Heat-space blowup for incomplete edge metrics 

The appropriate projective coordinates on ^1 are given as follows. Near the top 
corner of ff away from tf the projective coordinates are given by 

r ^ x ~ x y —y _ 

p = Vt, £ = -, £ = -, u = , y, z, z, 

P P P 
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where in these coordinates p, £, £ are the defining functions of the faces ff, rf and If, 
respectively. For the bottom corner of ff near If, the projective coordinates are given by 

t x y — y _ 

t = — , s = — , u = , x, y, z, z, 

X z X X 

where in these coordinates 7", Sj x circ the defining functions of tf, If and ff, respectively. 
For the bottom corner of ff near rf the projective coordinates are obtained by inter- 
changing the roles of x and x. The projective coordinates on near the top of td 
away from tf are given by 

T] = y/T, S = , U = -, Z = , x, y, z. 

7] 7] 7] 

In these coordinates tf is the face in the limit \(S, U, Z) \ — > oo, and ff and td are defined 
by x and rj, respectively. The blowup heat space is related to the original heat 
space M\ via the obvious 'blow-down map' (3 : ^j 2 — > M%, which in local coordinates 
is simply the coordinate change back to (t, (x, y, z), (x, y, z)). 

We can now state the asymptotic properties of the heat kernel as a polyhomogeneous 
distribution on the blowup These have been established by Mazzeo and the 

third author in [MV10]. The Friedrichs extension of the Laplace operator is a self- 
adjoint operator in the Hilbert space L 2 (M, vol {g)), and due to the specific form of the 
incomplete edge metric we have vol (g) = x^dx A vo\g dM (x) near the edge singularity. 
Hence one needs to take the singular factor of x^ of the volume into account in order 
to obtain the correct asymptotic behaviour. 

Theorem 2.2. [MV10] The heat kernel H lifts via the blowdown map P to a polyhomo- 
geneous distribution f3*H on with asymptotic expansion of leading order (— dim M) 
at the front face ff and at the diagonal face td, smooth at the left and right boundary 
faces, and vanishing to infinite order at tf . 

The heat kernel asymptotics were derived in [MV10] under a unitary rescaling trans- 
formation in the singular neighbourhood U = (0, 1) x dM C M of the edge singularity, 
similar to Briining-Seeley [BS88] 



$:C- t ([/)^C- t ([/), u^x" 2 



u. 



such that the transformed Laplacian foAo$ _1 is self-adjoint in L 2 {U, dx A vo\g dM (x)) . 
In particular the leading order of the heat kernel asymptotics in [MV10] refer to the 
heat operator as an integral operator with respect to the volume dx A vo\g dM (x) . This 
accounts for the difference in powers of the defining functions between Theorem 2.2 and 
[MV10]. 

To see this explicitly, we compare the actions of the heat operators with and without 
the unitary rescaling transformation. We may assume for simplicity that the heat kernel 
is compactly supported in U 2 x M + and find 



exp(-t$ o A o $ )u = / exp(-t$ o A o $ )(x,y, z,x,y,z)u(x,y,T)dxdvoldu(x) 

Ju 

x^ 2 exp(— tA)(x, y, z, x, y, z)x~^ 2 u(x, y, z)x^ dx dvobjM^) • 



u 
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Consequently we deduce the following relation between the heat kernels: 

exp(— 1$ o A o y, z, x,y,T) = (xx)^ 2 exp(— tA)(x, y, z, x,y,T). 

One may now check in projective blowup coordinates that the lifts of the heat kernels 
to the blowup heat space are related as follows 

(2.1.1) /T exp(-t$ o A o S- 1 ) = p£(p rfPlf // 2 /3*exp(-tA), 

where, e.g., p r f is a defining function for the right face. The rescaled heat kernel 
exp(— £$ o A o has been discussed in [MV10] and lifts to a density on the heat- 
space blowup with asymptotic expansion of leading order (—1 — 6) at ff, of leading order 
(//2) at If and rf, and of leading order (— dimM) at td. Employing the relation (2.1.1) 
we arrive at the actual statement of Theorem 2.2. 

Remark 2.3. A similar situation arises in the special case of a cone, where the edge 
B is collapsed to a point. This setup was considered by Mooers in [Moo96], where the 
heat kernel was considered under a unitary rescaling transformation as well. Without 
the rescaling one needs to take the singular factor of the volume x^dx A volg dM (x) 
into account, leading to different asymptotic behaviour. It appears this was not done 
in [JL03]. 

We will also need a more refined result on the asymptotic behaviour of the heat kernel 
at the right boundary face of M 2 . 

Proposition 2.4. Let (M m ,g) be a feasible edge space. Then the lift f3*H is a poly- 
homogeneous distribution on and letting s denote a defining function of rf, the 
expansion of f3*H takes the form 

f3*H ~ s°a (f3*H) + 0(s 7 °) as s 

for some 70 > ; where a (f3*H) is constant in the first fibre variable z. 

Proof. Recall the heat kernel construction in [MV10]. The initial approximative 
parametrix for the solution operator of C = (d t + A) is constructed in [MV10] by 
solving the heat equation to first order at the front face ff of ^ 2 . The restriction of 
the lift (3*(tC) to ff is called the normal operator Ng(tC) at the front face and feasi- 
bility of g, more precisely Definition 1.2 (i), ensures that N s (t£) is given in projective 
coordinates (r, s, u, z, x, y, z) by 

N s (tC) = r(d T - d 2 - fs-% + s- 2 A F , z + Af ) =: r(d T + Af<*> + A?). 

Here ^{F) = Kj~ xF z denotes the model cone. Since Ng(tC) does not involve derivatives 
with respect to (x,y,T), it acts tangentially to the fibres of the front face. Searching 
for an initial parametrix Hq, we solve the heat equation to first order at the front face, 
and want 

N s (tC o H ) = N s (tC) o N S (H ) = 0, 

which is the heat equation on the model edge < tf(F) x R*. Consequently, the initial 
parametrix H is defined by choosing Ng(H ) to equal the fundamental solution for the 
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heat operator Ns(tC), and extending Ns{H ) trivially to a neighbourhood of the front 
face. Using the projective coordinates (r, s, u, z, x, y, z) near ff we have 

(2.1.2) N S (H ) :=H v ^(t,s,z,1,T)H r \t,u,0), 

where H Rb denotes the euclidean heat kernel on M. b and is the heat kernel for the 

exact cone ^(F), as studied in [Che83], [Les97] and [Moo96]. Thus, the index set 
E C [0, oo) x {0} for the asymptotic behaviour of H Q at the right and left boundary 
faces is given in terms of the indicial roots 7 for the Laplacian Ag F ^ on the exact cone, 
i.e., 

H ~ s 7 a 7 (P ) as s — 7- 0, 

with the leading coefficient ao(H ) being harmonic on fibres and hence constant in z. 
Feasibility of g, more precisely Definition 1.2 (ii), ensures that the indicial roots are 
independent of the base point b G Y\ hence the index set E is discrete and H is indeed 
polyhomogeneous on The error of the initial parametrix Hq is given by 



P*{tC)H = \p*{tA g ) - rAj W - rAl 



Hn =: Pn 



The asymptotic expansion of H as s — > starts with s°a with no s°(log s) k , k > 1 
terms. The leading term s°a is annihilated by the edge vector fields sd s and 
d z since a is constant in z, and its order is raised by sd u . Consequently, since 

/3*(tA g ) — rAf^ — rAjf is of higher order in s, we find 



00 



P ~ s ' a '( P °) + sJ ' 1+la iAPo) as s -> 0, 

1=0 76_E* Z=0 

where E* := E\{0}. 

The next step in the construction of the heat kernel involves adding a kernel J to Hq, 
such that the new error term is vanishing to infinite order at rf. In order to eliminate 
the term s 7 a 7 in the asymptotic expansion of Pq at rf, we only need to solve 

(2.1.3) (-<9 S 2 - fs^ds + s- 2 A FjZ )u = s 7 (r" 1 a 7 ). 

This is because all other terms in the expansion of tC at rf lower the exponent in s by at 
most one, while the indicial part lowers the exponent by two. The variables (r, u, x, y, z) 
enter the equation only as parameters. The equation on the model cone is solved via 
the Mellin tranform and the solution u is polyhomogeneous in all variables including 
parameters, and is of leading order (7 + 2). Hence Hi = H + J expands near rf as 

00 00 
H x ~ s°a (#i) + s l ai{Hi) ^ +l a^(Hi) as s 0, 

1=2 7£E« 1=0 

where a (Hi) = a (H ) is constant in z. 

In the following correction steps the exact heat kernel is obtained from Hi first by 
improving the error near td and then by an iterative correction procedure, adding terms 
of the form if 1 o(P 1 ) fc , where Pi := tCHi is vanishing to infinite order at rf and td. Since 
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the leading coefficient a (Hi) is constant in z, we deduce d z Hi ~ 0(s 7 °) as s — > 0, where 
7o := min{2,7 e E*}. Consequently, we also have d z ii\ o [P\) k ~ 0(s 7 °) as s — >• 0. We 
find 

oo oo 

(3*H ~ s°a (f3*H) + s l ai((3*H) ^ +l a^i{P*H) as s ->• 0, 

where the leading coefficient ao((3*H) = ao(-ffo) is still constant in z. □ 

2.2. Heat kernel on spaces with a complete edge metric. We next describe 
the blowup space and asymptotic structure for the heat kernel of the Laplacian of a 
complete edge metric. The blowup proceeds exactly as in §2.1 except that in the first 
step we blow up the heat space M\ along the submanifold 

A = {(t, (0, y, z), (0, y, z)) E R + x («9M) 2 : y = y } C M 2 h . 

The final heat-space blowup for a complete edge metric is depicted in Figure 2. 




Figure 2. Heat-space blowup for complete edge metrics 



We now describe some coordinates on Near the left hand corner where the ff, 
tf and If meet we may use projective coordinates for valid away from x = 0: 

(2.2.1) r = Vt, s = — , u = - — — , x, y, z, z. 

x x 

Observe that in these coordinates, s and x are defining functions for If and ff, respec- 
tively. Near the top of the blown up diagonal, where td and ff meet, and away from tf 
we may use the coordinates 

(2.2.2) n = VI, S = U = V^y, Z = x, y, z. 

In these coordinates, rj is a defining function for td and x is a defining function for ff. By 
abuse of notation, we will denote the blowdown map with respect to any coordinates by 
(3, and will denote the composition of (5 with projections onto the left and right factors 
in Ml by (3 L : J^l -»■ K+ x M and (3 R : Jtl -> M, respectively. 

The asymptotic structure of the heat kernel of the Laplacian of a complete edge 
metric was described by Pierre Albin in [Alb07]. 
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Theorem 2.5. [Alb07] The heat kernel H lifts via the blowdown map (3 to a poly- 
homogeneous distribution (3*H on ^ 2 with asymptotic expansion smooth up to the 
boundary of the front face ff, with leading order (— dimM) at the temporal diagonal 
face td, and vanishing to infinite order at left, right and temporal boundary faces, If, rf, 
tfi respectively. 

Albin worked in the context of the heat operator acting on half-densities. We dispense 
with half-densities and instead regard the heat kernel as acting on functions (with 
compact support near the edge) by: 

tA/-/ „. _\ / tt u / „. _\ /- ~ ~\\ tt~ ~ ~\~-b-l. 



e f(x,y,z) = H (t, (x, y, z), (x, y, z)) f(x, y, z)x dxdvo\ 9M (x), 

J M 

where dvol^M denotes the volume with respect to the metric 4>*g B (x) + g F (x) on DM. 
This accounts for the difference in powers of the defining functions from [Alb07]. We 
sketch how to apply Albin's result to get the asymptotics of the heat kernel acting on 
functions instead of half-densities. 

We identify functions and half-densities on R + x M 2 , on the factor M. + x M, and on 
the factor M by 1 

f(t, (x, y, z), (x, y,z)) -H- f{t, (x, y, z), (x, y, ) x~ ( 2 ) (x) ~ ( 2 | dxdydzdxdydzdt | , 

_ (6+1) i /r, 

f(t,(x,y,z)) -H- f(t,(x,y,z))x ~^~\dxdydzdt\ 1 , 

_ (6+1) i /r, 

f(x,y,z) f(x,y,z)x 2 \dxdydz\ ' . 



As discussed in [Alb07], an element of the heat calculus A e \1> e ' Heat has an integral 

(6+/) (6+1) 

kernel that may be written as p td 2 p s 2 k ■ v, where pf ace is a defining function for 
the specified face; k is a function that vanishes to infinite order at If, rf, and tf, and is 
smooth up to the boundary at td and ff; and v is a smooth section of Vt 1 / 2 ^^^). Such 

2 

an operator A e \1> e ' Heat acts on half-densities by 



A f(t,(x,y,z))x ^ \dxdydzdt\ 1/2 

( (6+/) (6+1) (b+l) . .„ 

(2.2.3) =(Pl)Jp u 2 P ff 2 ku-{p R )*U\x,y,z){x)-— \dxdyd^ 2 ] 



To relate these half-densities, let us work in coordinates (2.2.2). We may take v 
\dSdUdZdxdydzdrj\ l l 2 . Pulling back our standard half-density on M 2 x M + , we find 



.(b+n (t±i) n/2,. 



f3*(x 2 (x) 2 \dxdydzdxdydzdt\ 

2 ("I) 

I— , s m , (6+1) (6+1) 6+/+2 

(2.2.4) =^(-1)^(1-^) — 2 x — —r}-^v. 



x-^(x(l - riS^-^i-l^^V^x^n^ldSdUdZdxdydzdnl 1 / 2 



1 Thcse densities have a factor of y'det g in appropriate places. In what follows we omit smooth 
factors. 
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Since 1 — r/S = - and we are working near the intersection of td and ff, the factor 
V / 2(— 1)^(1 — rjS)~^~2^ is smooth and uniformly bounded, so we absorb it into k. 
Denote a R := P R {(x) 2 \dxdydz\ 1/2 ). Using (2.2.3) and (2.2.4) we compute 



[A f](t, (x, y, z))x ^ \dxdydzdt\ 1/2 

( (t+J) (6+1) ^ ^ ^ (6+1) ^_ _ n /„ 

= {Pl)*[V 2 x 2 kv-(p R )*(f(x,y,z)(x) 2 \dxdydz\ 1/2 ) 



^ x {± ^Kv,(S,U,Z),(x,y,z)) 



(b+l) (b+f+2) (6+1) _ (6+1) ~ ~. ~ i /«j \ 

■x^~r] 2 /3*( x 3— (x) ~\dxdydzdxdydzdt\ ' )((3* R f)cr R J 

= (&). {v' m k(v, (S, U, Z), (x, y, *))(/%/) 

•/3*(x~ 2 (J) ~ 2 | dxdydzdxdyd'zdt \ 1 ^ 2 ) cr^ j 

= / # ( tV2 ' ^T/l' 7172^' x ' 2/> z ) /(*> ^ z)x- {b + X) dxdydz 

■ x'^—^xdydzdt^ 2 . 

Canceling the half density from the first and last components of this equation gives the 
action of the heat kernel on functions, with the stated asymptotics. 

2.3. Stochastic completeness of the heat kernel. Several of our arguments are 
simplified by the fact that the edge heat kernel is stochastically complete, i.e., satisfies 



(2.3.1) / H(t, (x,y,z), (x,j/,2))dvolg = 1, for all (x,y,z) e M,t > 0. 

Jm 

In the incomplete case this is a consequence of uniqueness of solutions to the heat 
equation. More precisely, with A being a self-adjoint unbounded operator on the Hilbert 
space L 2 (M, vol(g)) with domain @(A), the solutions to the initial value problem 

d t u + Am = 0, u(0) =u E @(A) C L 2 (M, vol(g)) 

are unique and in fact given by u(t) = e~ tA u E 3>(A) for any t > 0. Consequently, 
u = 1 is the unique solution in S>{A) to the heat equation with initial value 1, given in 
terms of the heat operator by u = 1 = e~* A l. This is precisely (2.3.1). 

In the complete edge case, u = 1 does not lie in the Hilbert space L 2 (M, vol (g)) 
and on general complete manifolds, solutions of the heat equation need not be unique. 
However by a result of Yau [Yau78] we have (2.3.1) on complete manifolds with Ricci 
curvature bounded from below; see [Hsu89] for a discussion of this and related results. 
In our setting, since an edge metric can be viewed at highest order as the product of an 
asymptotically hyperbolic metric and a compact metric, it is straightforward to check 
that the Ricci curvature is bounded from below. 

3. Mapping properties of the heat kernel on edge manifolds 

3.1. Estimates on spaces with incomplete edge metrics. We prove Theorem 1.5. 
The main task in the proof is to check the estimates, as the continuity of e~ tA f and 
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its derivatives are easy to obtain by the dominated convergence theorem in each of the 
coordinate systems that will be discussed. The estimates in Theorem 1.5 are classical 
for the lift of the heat operator to M. 2 h supported away from the front face. So we 
can assume without loss of generality that the lift of the heat kernel to the heat space 
blowup ^ 2 is compactly supported in an open neighbourhood of the front face. Then 
the heat operator acts on functions as in (2.0.1) and the proof amounts to estimating 
{k + n) edge derivatives V^ +n locally: 

(3.1.1) / V* +n H (t,(x,y,z),(x,y,z))f(x,y,z}x f dxdvol dM (x}. 

J M 

The estimates require a precise understanding of the asymptotic behaviour of H . Since 
the asymptotics of the heat kernel H near the front face are non-uniform, one needs to 
estimate the integral (3.1.1) for the integrand supported near the various corners and 
boundary faces of the heat space blowup We thus separate H into four compo- 

nents, with their lifts to the heat space blowup ^1 compactly supported near each of 
the three corners of the front face and near the diagonal of This corresponds to 
different asymptotic regimes as the variables (t, x, x) approach zero, with the associ- 
ated projective coordinates from §2.1 capturing the various asymptotic behaviour. We 
establish the estimates for the integral operators defined by each of the components 
separately. 

Estimates near the lower left corner of the front face 



We assume that the heat kernel H is compactly supported near the lower left corner 
of the front face. Its asymptotic behaviour is appropriately described in the projective 
coordinates 

t x y — y _ 
r = — , s = ~, u = , x, y, z, z, 

X z X X 

where r, s, and x are the defining functions of tf, If and ff, respectively. The coordinates 
are valid whenever (r, s, u) are bounded as (t, x, x) approach zero. The edge vector fields 
obey the following transformation rules: 

0*{xd x ) = -2rd T - sd s - ud u + xd x , (3*(xd y ) = d u + xd y , /3*{d z ) = d z . 

Note that we abuse notation by writing xd x , xd y vector fields on the right side of these 
expressions. By Theorem 2.2 we infer that for any n G No, 

P*(y?H)(T,x,y,z,s,u,z) = x~ m G{T,x,y,z,s,u,z), 

where G is bounded in its entries and infinitely vanishing at r goes to zero, with the 
same compact support as (3*H. For the transformation rule of the volume form we 
compute 

(3*(x^dx dvobjM^)) = h ■ s^x m ds dudz, 
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where h is a bounded distribution on Hence we can write for any / e C]?(M,g) 
V^e~ tA f{x, y,z)= / V^H (t, {x, y, z), (x, y, z)) f(x, y, z)x f dx dvo\ dM (x) 

J M 

= J h • G • f(sx,y — ux,z) ds dudz. 

Pulling out the supremum of / and estimating h ■ G against a constant, we find 

\V e n e- tA f(x, y,z)\ < C ll/IU < C-\\f\\ cb 

Note that this estimate also accounts for the second statement of the theorem, since 
C«(M,g)cC° e (M,g). 

Estimates near the lower right corner of the front face. 

Applying our argument for the left corner of the front face near the right corner 
appears to produce singular behaviour in powers of the defining function for the front 
face. However, the argument is not symmetric in the (x, y, z) versus (x, y, T) variables 
as these variables play different roles in the integration (2.0.1). Thus we may use a 
second change of variables to take advantage of the infinite order decay at the temporal 
face. To begin, we use the projective coordinates 

~ t _ x „ y — y 

T ^7, s — , u — , 2, x, y, z, 

X z X X 

where r, s, and x are the defining functions of tf, rf and ff, respectively. The coordinates 
are valid whenever (r, s, u) are bounded as (t, x, x) approach zero. The edge vector fields 
obey the following transformation rules: 

f3*(xd x ) = sdg, (3*{xd y ) = ads, f3*{d z ) = d z . 

By Theorem 2.2 we infer that for any n e N , 

P*(V?H)(t, s, u, z, x,y,z) = x- m G(r, s, u, z, x,y,z), 

where G is bounded in its entries and vanishing to infinite order in r as r goes to 
zero, with the same compact support as (3*H. Of course G depends on the particular 
derivatives taken, but we will not express this in the notation. 

In order to estimate (3.1.1) we again perform a change of variable, this time using 
the time variable in order to take advantage of the decay of the heat kernel at tf. This 
gives 

x^ dx dvohjM^) = h ■ T~ 1 x 1+ f +b dr du dz = h ■ r~ 1 x m dr dudz, 
where h is a bounded distribution on Therefore for any / e C^(M,g) 

V e n e- tA f(x,y,z) 

= J x~ m ([3*G) (t, (x,y, z), (x,y,z)) f(s/i/Vr,y — uVt/V^,z^jh ■ T~ 1 x m drdudz 
= J r^/^G) (t, (x,y, z), (x,y,z)) f{Vt/V?,y — uVt/V^,z)h ■ drdudz, 



16 ERIC BAHUAUD, EMILY B. DRYDEN, AND BORIS VERTMAN 

where G is vanishing to infinite order in r and bounded in its other entries. We conclude 

\\V:e- tA f{xiy,z)\\c«<C\\f\\ c o. 
Further, we find from the expressions for the edge vector fields that 

||e-* A /(x )2/ ,z)|| c ,<L7||/|| c o<L7||/|| Cefc 

and 

\\e-^f{x,y,z)\\ c ,<C\\f\\ c0 <C\\f\\ c? . 



Estimates near the top corner of the front face 

The asymptotic behaviour of the heat kernel H, compactly supported near the top 
corner of ff, is appropriately described in the projective coordinates 

r ^ x ~ x y — y ~ 
p = Vt, 4 = -, 4 = -, u = , y, z, z, 

P P P 

where p, £, and £ are the defining functions of the faces ff, rf and If, respectively. The 
coordinates are valid whenever (p, £, £) are bounded as (t, x, x) approach zero. The edge 
vector fields obey the following transformation rules: 

p*(xd x ) = 4% (3*{xd y ) = £d u + fidy, f3*{d z ) = d z . 

Hence by Theorem 2.2 we infer that for any n e N , 

/?* (V e n H) (p,Z,y,z,Zu,z)= p- m G(p, S,y,z,Zu,z), 

where G is bounded in its entries and has the same compact support as (3*H. For the 
transformation rule of the volume form we compute 

0*(x f dxdvol dM (x)) = h ■ I s p 1+}+b dldudz, 

where h is a bounded distribution on jM'^- Hence we can write for any / G C^(M,g) 



V™e tA f(x,y,z) = / V™H (t, (x,y,z),(x,y,z)) f(x,y,z)x f dxdvolgu(x) 

J M 

= J h-G-f(p£,y — up,z)^d£dudz. 

Pulling out the supremum of / and estimating h ■ G ■ against a constant, we find 

\V e n e- tA f(x,y,z)\< C||/|U< C7.ll/llc*. 
This estimate again accounts for the second statement of the theorem as well. 

Estimates where the diagonal meets the front face 

The asymptotic behaviour of the heat kernel with compact support in the neighbour- 
hood where td meets ff is appropriately described in the projective coordinates 

(3.1.2) rj - -, S - - ?r , U--^-,Z- —J*-, x, y, z 
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where in these coordinates 77 and x are the defining functions of td and ff, respectively. 
The edge vector fields obey the following transformation rules: 

F{xd x ) = -r)d n + -d s + Zd z + xd x , 0*{xd y ) = -d v + xd y , /3*(d z ) = -d z + d z . 

r] rj V 

Hence by Theorem 2.2 we infer that for any n e No, 

P*(V e n H)( V , S, U, Z, x, y, z) = x- m V ~ m - n G( V , S, U, Z, x, y, z), 

where G is bounded in its entries (in fact, infinitely vanishing as 1(5, U, Z)\ — > oo) and 
has the same compact support as ft*H. 

The transformation rule of the volume form is 

(3*(x f dx&vo\ aM (x)) = h ■ {xr]) m (l - r]S) f dSdUdZ, 

where h is a bounded distribution on Hence for any / e C®(M,g), we have 

V™e~ tA f(x, y,z) = / V™H (t, (x, y, z), (x, y, z)) fix, y, z)x f dx dvol 9M (x) 

J M 

= J T]- n (l - r)S) f h ■ G f(x(l - r]S),y - xrjU, z - rjZ) dSdUdZ. 

Pulling out the supremum of / and estimating x n (l — r]S)^h ■ G against a constant, we 
find 

\V:e- tA f(x,y,z)\<C-t- n / 2 \\f\\ c o. 

This completes the proof of estimate (1.0.1) for k = 0. 

For the general case, let / G C*(M, g) and consider the impact of k + n edge deriva- 
tives. Transforming the edge derivatives into projective coordinates near where the 
diagonal meets the front face, we see that the only possible singular behaviour comes 
from components of the form rj~ l ds,rj~ l du and r]~ 1 dz- We consider rj~ 1 ds (the others 
are similar) and estimate the integral 

F : — j V 1 (d s (3*H) (rj, S, U, Z, x, y, z)f ■ h(xr]) m (l - r]S) f dS dU dZ 

= J V~ X (dsG) (77, S,U, Z,x,y, z)f ■ h(l - rjSy 'dS dU dZ. 

We perform integration by parts, where the boundary terms lie away from the diagonal 
and hence are infinitely vanishing for t — > by the asymptotic behaviour of the heat 
kernel. Omitting these irrelevant terms we obtain 

F = J G{r), S,U, Z,x,y, z)r 1 - 1 d s (f ■ h(l - r]S) f ) dS dU dZ 

= J G[r]- l d s f {x{l-r}S),y-xr]U,z-r)Z)] h(l - r]S) f dS dU dZ 
+ j Gf- [q^dshixil - r]S),y - xr]U, z - rjZ) (1 - r]S) f ] dSdUdZ. 
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The second integral yields no ^ _1 -singular behaviour and is estimated as above for 
(k,n) = (0,0). For the first integral note 

r}~ l d s f -r)S),y- xrjU, z - r]Z) = - 1 (xd$f) (x(l -r]S),y- xrjU, z - r]Z) . 

(1 - 7]S) 

Since G(r], S, U, Z, x, y, z) is bounded in its entries, and in fact infinitely vanishing as 
\(S, U, Z) \ — > oo, we can estimate 

\F\<C\\f\\ ci . 

Repeating the same argument for k > 1 edge derivatives we find that the heat 
operator is bounded in its action on C^(M,g). In view of the proof for k = 0, this 
proves the estimate (1.0.1) in full generality and it remains to check (1.0.2). Recall 
from §2.3 that the heat operator satisfies e~ tA l = 1 by uniqueness of solutions to the 
heat equation and hence V 2 e~ tA l = 0. In particular, 

Ve e ~ tA f = / K 2 #(^ V-> z )i (^ y> z)x S dxfao\ dyi (x) 

J M 

V 2 H(t, (x,y,z),(x,y,z))[f(x,y,z) - f(x,y,z)]x f dxdvol dM (x) =: E. 

Using the same approach as above, we obtain 

E= I rf 2 G [f(x(l-r]S),y-xr)U,z-r)Z)- f(x, y,z)]h(l- r]S) f dS dUdZ. 
Jm 

Estimating the difference in terms of the Holder norm of /, we find 
\E\ < \\f\\c S J V~ 2 G- [d ((x - xr]S,y - xrjU, z - rjZ) , (x,y, z))] a h(l - riS) f dS dU dZ 

< C\\f\\c? J ri~ 2 Hl ~ vS) f G ■ (\xr]S\ a + \xr]U\ a + |(2 - r)S)xr]Z\ a ) dS dU dZ 

< C'WfWc? j V~ 2+a h(l - vS) f G( V , S, U, Z, x, y, z) (\S\ a + \V\ a + \Z\ a ) dS dU dZ. 

Since G(rj, S,U, Z,x,y, z) is bounded in its entries, and in fact infinitely vanishing as 
\(S, U, Z) \ — > oo, we can estimate 



\e\ < cr l+a/2 



We have proved Theorem 1.5 for the heat operator acting on C*(M,g). In order 
to complete the proof of Theorem 1.5 it suffices to prove that the image of C^(M,g) 
under the heat operator lies in ^ e fc+2 (M, g). This is a direct consequence of Propo- 
sition 2.4, which asserts ^-independence of the leading coefficient in the right face 
expansion of H. Indeed, for any fixed t > and any / e C^(M,g), Proposition 
2.4 implies {d z {V x ■ ■ ■ V k+2 )e- tA f){0, y, z) = for any Vj G V e ,j = l,..,k + 2. Hence 
e -tA (pl(M,g)) C ^ +2 (M,g). This completes the proof. 
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3.2. Estimates on spaces with a complete edge metric. In this section we prove 
Theorem 1.6. Note that in contrast to the incomplete edge case, the heat operator 
here respects the additional weighting of the edge spaces. Consequently we derive more 
refined estimates. 

Once again the main task in the proof is to check the estimates, as the continuity 
of e~ tA f and its derivatives are easy to obtain by the dominated convergence theorem 
in each of the coordinate systems that will be discussed. The proof of the mapping 
properties amounts to estimating edge derivatives of the heat kernel lifted to these 
estimates are made locally in various coordinate patches of 

The estimates for the heat kernel supported near the diagonal away from the front 
face or supported in the interior of jft"^ are known (e.g., [LSU67]), so our discussion 
reduces to the case of the heat kernel compactly supported in an open neighbourhood of 
the front face. Within this neighbourhood, we consider the heat kernel being supported 
in each of three coordinate charts: one near each of the left and right corners of the front 
face, for which we use coordinates (2.2.1), and one near the intersection of td and ff, for 
which we use coordinates (2.2.2). Although the argument is not completely symmetric 
in the left and right corners due to the different roles played in the integration by tilded 
and untilded variables, the estimation is analogous, as we indicate below. These three 
charts correspond to different asymptotic regimes as the variables (t, x, x) approach 
zero. 

Before we specialize to the coordinate charts, we note that the estimates of Theorem 
1.6 for weighted spaces follow from the estimates for unweighted spaces. Indeed, proving 
the estimates for weighted spaces is equivalent to proving the unweighted estimates for 
the conjugated operator X~ w e~ tA X w , where X refers to the obvious multiplication 
operator. Written as an integral operator, we have 

y,,)=[ X -H it, (,, (x, y,Ii) i» • /ft y, ^-VJdvoMJ). 

J M 

One can check that in the coordinate systems (2.2.1) and (2.2.2), the conjugated integral 
kernel x~ w Hx w has the same asymptotics as H at all boundary hypersurfaces of 
Thus we may assume w = without loss of generality. 

Estimates near the lower left corner of the front face 

We begin by assuming that the heat kernel is supported near the left-hand corner 
and work in the coordinates (2.2.1) 

r x y — y _ 
r = vt, s = — , u = , x, y, z, z. 

x x 

In these coordinates r, s, and x are the defining functions of tf, If and ff respectively. 
The coordinates are valid whenever (r, s, x) are bounded as (t, x, x) approach zero. The 
edge vector fields obey the following transformation rules: 

0*(xd x ) = -sd s - ud u + xd x , f3*(xd y ) = d u + xd y , /3*{d z ) = d z . 

Hence by Theorem 2.5 we find that for any n G No, (3*(V™H)(t, s, u, x, y, z, z) is bounded 
in x and vanishing to infinite order as r — > and s — > 0. For the transformation rule 
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of the volume form we compute 

f3* (x~ b ~ l dx &vo\qm{x)) = h ■ s~ b ~ l ds dudz, 
where h is a bounded distribution on Hence for any / e C^(M,g), we have 

Pulling out the supremum of /, we find 

IK^/foy,*)! < WfWooJ h-P*{V?H)s- b - l dsdudz. 

Since /3*(V^ n i7)(r, s, w, x, y, z, 2) is vanishing to infinite order as s — > 0, we can estimate 
the integral against a constant and find 

\V:e-^f{x,y,z)\<C\\fU<C .\\f\\ cb 

for any k e No- The second statement of Theorem 1.6 follows as well, since C°(M, g) C 
C° e (M,g). 

Estimates near the lower right corner of the front face 

In the complete case the estimates near the right corner are simplified by the infinite 
order vanishing of the heat kernel at the right face. We begin by assuming that the 
heat kernel is supported near the right-hand corner and work in the coordinates (2.2.1) 

r ~ x ~ y — y ~ ~ ~ 

t = vt, s = —, u — — — — , x, y, z, z. 

In these coordinates r, s, and x are the defining functions of tf, rf and ff respectively. 
The coordinates are valid whenever (r, ~s, x) are bounded as (t, x, x) approach zero. The 
edge vector fields obey the following transformation rules: 

(3*(xd x ) = sc\, (3*(xd y ) = sdu, 0*(d z ) = d z . 

Hence by Theorem 2.5 we find that for any n G No, f3*(V™H){r, s, u, x, y, z, z) is bounded 
in x and vanishing to infinite order as r — > and If — > 0. For the transformation rule 
of the volume form we compute 

/3*{x~ ~ x dx dvolaM(a^)) = h ■ ds du dz ', 
where h is a bounded distribution on M. 2 h . Hence for any / e C*(M,g), we have 



V™e tA f(x,y,z) = J h-P*(V™H)(T,s,u,x,y,z,J)f(x/s,y — xu,z)s 1 dsdudz. 

Pulling out the supremum of /, we find 

\V,?e- tA f(x,y,z)\ < \\fU J h- ^{V^H^dsdudz. 

Since (3*(V™H)(t, s, u, x, y, z, z) is vanishing to infinite order as s — > 0, we can estimate 
the integral against a constant and find 

IK n ^ A /(^y^)l<^ll/lloc<G-||/|U, 
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for any k e No- The second statement of Theorem 1.6 follows as well, since C°(M, g) C 
C° e (M,g). 

Estimates near the diagonal 

The asymptotic behaviour of the heat kernel with compact support in the neigh- 
bourhood where the diagonal meets the front face is appropriately described in the 
projective coordinates (2.2.2) 

n Q x-x y-y z-z 

<q = Vt, S = ■=-, U = 7T , Z = — — , x, y, z. 

xyt xyt yr 

Recall that in these coordinates r\ and x are the defining functions of td and ff, respec- 
tively. The edge vector fields obey the following transformation rules: 

(3*(xd x ) = (1 - V S)-d s - Ud v + xd x , 0*(xd y ) = -d v + xd y , f3*(d z ) = ~d z + d z . 

7] 7] T] 

Hence by Theorem 2.5 we infer that for any n e N , 

P*(V e n H)(r), S, U, Z, x, y, z) = ^^G^, S, U, Z, x, y, z), 

where G is bounded in its entries, and in fact infinitely vanishing as \(S, U, Z)\ — > oo, 
with the same compact support as (3*H. The coordinates (x,y,z) on the second copy 
of M are expressed in terms of the projective coordinates (2.2.2) by 

x = x(l — rjS), y = y — xrjU, z = z — rjZ. 

Thus the transformation rule of the volume form is 

^{x-^alxdvolduix)) = hrPil-riSy^dSdUdZ, 

where h is a bounded distribution on Hence for any / 6 C®(M,g), we have 

D : = / V* l H (t, (x, y, z), (x, y, z)) f(x, y, z)x~ b ~ l dx &vo\ dM (x) 
Jm 

r)- n {l-r)Sy b - l h-Gftx{l-T)S),y-xr)U,z-r)Z) dSdUdZ. 

Pulling out the supremum of / and estimating (1 — r]S)~ b ~ 1 h ■ G against a constant, we 
find 

D<c-r n / 2 \\f\\ c o. 

This completes the proof of estimate (1.0.3) for k = 0. For the general case, let 
/ G C*(M,g) and consider the impact of (k + n) edge derivatives. Transforming the 
edge derivatives into projective coordinates near the intersection of the diagonal and 
the front face, we see that the only possible singular behaviour comes from components 
of the form r/ -1 <9s, n~ l djj and n~ l dz- We consider rj~ 1 ds (the others are similar) and 
estimate the integral 

F : = J r,- 1 {d s f3*H) (rj, S, U, Z, x, y, z)f ■ hrf\\ - r/Sy^dS dU dZ 
r]- 1 (d s G) (v, S, U, Z, x, y, z)f ■ h{l - rjSy^dS dU dZ. 
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We integrate by parts and note that G is infinitely vanishing for S — >■ ±00 by the 
asymptotic behaviour of the heat kernel. Omitting this boundary term we obtain 

F = J G{r],S,U,Z,x,y,z)7]- 1 d s {f ■ h(l - r ] S)- b - 1 ) dS dU dZ 

G [v~%f (z(l -r}S),y~ %VU, z - r]Z)] h{l - rjSy^dS dU dZ 
+ j Gf- [T]- l d s h(x{l -r)S),y-xr)U,z-r]Z) (1 - r/S)^ 1 ] dSdUdZ. 

After differentiating in the second integral we see that no r]~ 1 term remains, and this 
integral can be estimated using the same approach as for D above. For the first integral 
note that 

V^dsf (x(l -r)S),y- xrjU, z-r)Z) = - 1 (xd s f) (x(l -r]S),y- xrjU, z - r]Z) . 

(1 - T}S) 

The edge derivative xdxf is bounded above by H/Hc 1 ! the remaining terms in the 
integral can be estimated as for D since G(rj, S, U, Z, x, y, z) is bounded in its entries 
and in fact infinitely vanishing as \(S, U, Z) \ — > 00. Thus |F| < C||/||ci- 

Repeating the same argument for k > 1 edge derivatives we find that the heat 
operator is bounded in its action on C^(M,g). In view of the proof of estimate (1.0.3) 
for k = 0, the first statement of Theorem 1.6 holds in full generality and it remains 
to check the second statement. Recall from Section 2.3 that the heat operator satisfies 



-tA-i 



V e 2 e~ tA f 



1 by stochastic completeness, and hence V. e 



0. In particular, 



V e H(t, (x,y,z), (x,y,z))f(x,y,z)x dx dvo\ 9M (x) 

V?H(t, (x,y,z),(x,y,z)) [f(x,y,z) - f(x, y, z)\ x~ b ~ l dx dvol 9M (^) =: E. 



M 



Using the same approach as for D above, we obtain 



E 



r]- 2 G [f (x(l -r}S),y- xr]U, z-rjZ)- f(x, y, z)} h(l - rjSy^dS dU dZ. 
Estimating the difference in terms of the Holder norm of /, we find 
\E\ < ||/|| c - I r]- 2 G- d((x - xr]S,y - xr]U, z - rjZ) ,(x,y,z)) a h(l - rjSy^dS dU dZ 



<c\\f\\cs v H^-vsr^G. 



T]S 



(2 - V S) 



T]U 



T)S) 



\r)Z\ a dSdUdZ 



c? 



-2+a 



G(r),S,U,Z,x,y,z)(\S\ a + \U\ a + \Z\ a )dSdUdZ. 



Since G(r), S,U, Z,x,y, z) is bounded in its entries, and in fact infinitely vanishing as 
\(S, U, Z) \ — > 00, we can estimate 

\E\<Ct- i+a / 2 \\f\\ c? . 
The proof of Theorem 1.6 is now complete. 
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4. Short- time existence of solutions to semilinear parabolic equations 

In this section we provide an application of our estimates to prove short-time existence 
of solutions to some semilinear parabolic equations on both complete and incomplete 
edge spaces. We follow [JL03, JL04] closely. The underlying idea is based on [Tay96]. 

Consider an equation of the form 

dfU — Au = Q(u,V e u), 
w(0) = M , 

where A denotes the appropriate edge Laplacian and Q denotes inhomogeneous terms 
that involve up to one edge derivative of u. Applying Duhamel's principle, we transform 
the problem to an equivalent integral equation: 



(4.0.1) 



u = e- t/x u + / e-V—>*Q(u, V e u)ds. 
Jo 

With appropriate assumptions and the estimates of Theorems 1.5 or 1.6, this equation 
may be solved by the contraction mapping principle, as we now explain. 

Recall that a semigroup S t of operators on a Banach space X is strongly continuous 
if for all / G X, 

llSJ-Zlbr— >0ast— >0. 

Taylor proves: 

Theorem 4.1. [Tay96, Proposition 15.1.1] Suppose that X andY are Banach spaces 
such that: 

(i) e~ tA : X — > X is a strongly continuous semigroup, for t > 0. 

(ii) Q : X — > Y is locally Lipschitz. 
(hi) e~ tA : Y — ► X, fort > 0. 

(iv) For some 7 < 1, 

||e- tA |U ( y,x)<cr^. 

Then for any uq G Y , the initial value problem (4.0.1) has a unique solution u G 
C([0, T], X), for some T > 0. T may be estimated from below in terms of \ \f\\x- 

In view of this theorem and our previous estimates, it remains to prove strong conti- 
nuity of the semigroup e~' A on either X = ^ e fc (M, g) or X — C*(M, g), as appropriate, 
to deduce short-time existence of solutions to (4.0.1) for certain locally Lipschitz Q. 
See [Gy09] for the classical proof. 

Proposition 4.2. Let (M, g) be an m- dimensional Riemannian manifold with a feasible 
incomplete non-iterated edge metric. Let e~* A denote the heat operator corresponding to 
the Friedrichs extension A of the associated Laplacian on (M, g). Then e" <A is strongly 
continuous on ^(M,g), for any k > 0. 

Proof. We prove the statement by adapting the classical proof of strong continuity of 
the heat operator on closed (non-singular) manifolds to the present setup. Assume first 
k = 0. Using stochastic completeness of the heat kernel we find 

(e- tA f)(p,t)-f(p)= [ H(t,p,p)(f(p)-f(p))dvo\ g (p). 

JM 
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Note that / e ^(M,g) and hence for any e > there exists some 5(e) > 0, such 
that for d(p,p) < 5(e) one has \f(p) — f(p)\ < e. For any given e > we separate the 
integration region into 

M+ := {p\d(p,p)>5(e)}, 
(4.0.2) M~ := {p\ d(p,p) < 5(e)}. 

Employing continuity of / we find 

\e~ tA f ~f\ = \J M H(t,p,p) (f(p) - f(p)) dvol g (p)| 

< j M+ H (t,p,p) \f(p) - f(p)\ dvol g (p) + f M _ H (t,p,p) \f(p) - f(p)\ dvo\ s (p) 

< 2^J H/llo f M+ H (t,p,p) « dvo\ g (p) + e J M _ H (t,p,p) dvol g (p). 

The second integral is bounded by e. We must verify in projective coordinates that 
the first integral is bounded, uniformly in (t,p, e). It suffices to assume the heat kernel 
is supported near the front face, since away from ff the estimates reduce to the clas- 
sical case of a manifold without boundary near the diagonal. Assume first that the 
heat kernel is supported in a neighbourhood where td meets ff, and use the projective 
coordinates (3.1.2): 

In these coordinates rj and x are the defining functions of td and ff, respectively. By 
Theorem 2.2 we infer 

P*H( V , S, U, Z, x, y, z) = x- m V ~ m G( V , S, U, Z, x, y, z), 

where G is bounded in its entries, and in fact infinitely vanishing as 1(5", U, Z)\ — > oo, 
with the same compact support as j3*H. The transformation rule for the volume form 
is 

f3*(x f dxdvo\ 9M (x)) = h(xr]) m (l - 7]S) f dSdUdZ, 

where h — h(r), x(l — r/S),y — xr/U, z — r/Z, x, y, z) is a bounded distribution on ^M^. 
We may estimate d(p,p)/\/t: 

d(p,P ) 
Hence we find 

/ H(t,p,p)^jB-dyol g (p) <C [ h(l - v S) f G ■ V\S\ 2 + \U\ 2 + \Z\ 2 dSdUdZ < C", 
J m+ yt J 

where C" is independent of (t, x, y, z) and e. Entirely analogous estimation works near 
the top corner of the front face. In the estimate near the left corner (resp. right corner) 
of the front face one finds that d(p,p)/\/i is bounded by a power of r _1 (resp. r _1 ), 
which may be absorbed into the heat kernel as it vanishes to infinite order as r — > 
(resp. r — > 0) at this corner. Therefore we obtain 

\\e- tA f-f\\o<C"^\\f\\ + e. 



< cV|s| 2 + |c/| 2 + |z| 
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Thus, for any given e > we can estimate \\e tA f — f\\ < 2e for \ft < e5(e)/(C" 
This proves strong continuity of the heat operator on ^(M,g). 

We prove strong continuity on %f e (M, g>) for k > 1 in a similar fashion. Near the 
diagonal it requires an integration by parts argument like the one in §3.1. The esti- 
mates away from the diagonal are the same as for ^(M,g), so we assume that (3*H 
is compactly supported in a neighbourhood where td meets ff. The edge vector fields 
obey the following transformation rules: 



1 11 

(xd x ) = -r)d„ + -d s + Zd z + xd x , (3*(xd y ) = ~d v + xd y , /3*(d z ) = -d z + d z . 

7] 7] 7] 



We consider [^^(e tA f — /)||o- Using stochastic completeness of the heat kernel, we 
find 



F := xd x (e tA f - f)= J (xd x H)f(x,y,z)x f dxdvo\ 
- J (xd x )[Hf(x,y,z)x f dxdvo\ dM {x)] =: F 1 - F 2 



dM{X) 



Next we transform to projective coordinates and integrate by parts in S, where the 
boundary terms lie away from the diagonal and hence are infinitely vanishing for t — > 
by the asymptotic behaviour of the heat kernel. Omitting these irrelevant terms, we 
obtain 



F x = J \- n d T , + ^d s + Zd z + xdA [(xr))- m G(r),S,U,Z,x,y, 
x f (x(l - 7]S),y - X7]U, z ~ r]Z) h(xr]) m (l - 7]S) f dS dU dZ 
= J [(-vdr, + Zdz + xd x )(xri)- m G]- fh(xri) m (l-riS) f dSdUdZ 

G \(-ds) fl h(l - rjSy 'dS dU dZ 

(xrj)- m G-f ' ' 



-d s ) h{xr]) m {l-7]S) f 



dS dU dZ. 
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We perform similar computations for F^. 



(xd x H)f(x,y,z) + H(xd x f) x f dx dvol dM (x) 



-r]d v + -d s + Zd z + xd x 
V 



{xrf)- m G) f ■ h(xri) m (l-riS) f dSdUdZ 



+ J G(r],S,U,Z,x,y,z)(xd x f(x,y,z))h(l-r]S) f dSdUdZ 

= J [(-vd ri + Zd z + xd x )(xri)- m G]- fh(xri) m (l-riS) f dSdUdZ 

- J (xr])- m G ■ f Qd s ^ h(xr]) m (l - r]S) f dS dU dZ 

+ J G(ri,S,U,Z,x,y,z)(xd x f(x,y,z))h(l-riS) f dSdUdZ. 

Thus F = F 1 — F 2 becomes 

F = J [(-r ] d v + Zd z + xd x )(xr ] )- m G(r ] ,S,U,Z,x,y,z)jh(xr ] ) m (l-r ] Sy 
x (f(x(l — r)S),y — xr/U, z — r/Z) — f(x, y, z))dS dU dZ 



I 



G(r), S, U, Z, x, y, z) 



-d s )h-(i- v sy 

T] 



x (f(x(l — r)S),y — xrjU, z — r]Z) — f(x, y, z))dS dU dZ 



G 



-d s f(x(l -r]S),y- xrjU, z - rjZ) + xd x f(x, y, z) 



h{l-r]S) f dSdUdZ. 



Now, each of the three integrals is estimated as above for k = by separating the 
integration region into M+ and M~ for any given e > 0. Note that in the final integral 
we use the fact that / G ^(M, g). Higher order and other edge derivatives may be 
estimated in a similar way. □ 

Proposition 4.3. Let (M, g) be an m- dimensional Riemannian manifold with a feasible 
complete edge metric. Let e~ tA denote the heat operator of the unique self-adjoint 
extension A of the associated Laplacian on (M,g). Then e~ tA is strongly continuous 
on C*(M,g), for any k > 0. 

Proof. We proceed as in the incomplete edge case, again modeling the proof after the 
classical estimate on a closed manifold. We estimate near td and ff; the estimation near 
either the left corner or right corner is simpler than what follows due to the infinite 
order vanishing of the heat kernel at tf. 

We first consider the case k = and show \\e~ tA f — /|| — > as t — > 0. Using 
stochastic completeness of the heat kernel we find 

(e- tA f)(p,t)-f(p)= [ H(t,p,p)(f(p)-f(p))dvol g (p). 

J M 
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Recall that / G C® (M) means / is continuous on a compact manifold with boundary, 
so for every e > we obtain a uniform 6 = 5(e), where d(p,p) < 5 implies \f(p) — f(p)\ < 
e. We separate the integration region into M e + and M~ as in (4.0.2) and find 



-tA 



/-/I 



H(t,p,p)(f(p)-f(p))dvo\ g ( 

M 



< / H(t,p,p)\f(p)-f(p)\dvo\ g (p)+[ H(t,p,p)\f(p)-f(p)\dvo\ 

J M+ J M~ 

< 2^||/||o / ^(t,P,p)^dvol g (p) + e . 

°l e J Jm+ yt 

Assume that the heat kernel is supported in a neighbourhood where td meets ff, and 
use the appropriate projective coordinates (2.2.2) 

n Q x-x y-y z-z 

rj = Vt, S = ■=-, U = 7T , Z = — —, x, y, z. 

xyt xyt yt 

In these coordinates rj and x are the defining functions of td and ff, respectively. The 
edge vector fields obey the transformation rules: 

f3*(xd x ) = (1 - V S)-d s - Ud v + xd x , P*(xd y ) = -d v + xd y , f3*{d z ) = ~d z + d z . 

7] 7] 7] 

Hence, by Theorem 2.5 we infer 

/3*H(t], S, U, Z, x, y, z) = r]~ m G(r], S, U, Z, x, y, z), 

where G is bounded in its entries, and in fact infinitely vanishing as \(S, U, Z)\ — > oo, 
with the same compact support as f3*H. The coordinates (x, y, z) on the second copy 
of M are expressed in terms of the new projective coordinates by 

x = x(l — r/S), y = y — xr/U, z = z — r\Z. 

From there we compute the transformation rule of the volume form 

f3*(x- b - 1 dxdvol dM (x)) = hr] m {l-r]S)- h - l dSdUdZ, 

where h = h(r),x(l — r]S),y — xr]U,z — r]Z,x,y,z) is a bounded distribution on 
whose arguments we will suppress. Further, we can estimate 



^<J ro S \, 2 + (0 u \, 2 + z*<cvsi + ui + zi. 

v/t ]j(2-7]S) 2 (2-r)S) 2 
Hence we can write after cancellations 

|e- tA /-/|< 2-^||/||o / H(t,p,p)^Pldvo\ g (p) + e 
5(e) Jm+ Vt 

< CtttII/IIo f HI - vS)-*' 1 -G-VS 2 + U 2 + ZHS dU dZ + e. 
^(e) J 

<^||/||o + £ , 
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where C is independent of (p,t) (see §3.2 for similar estimation near where td meets 
ff). As in Proposition 4.2, we conclude 

lim||e- iA /-/||o = 0, 

completing the case k = 0. 

To prove strong continuity on C^(M,g) for k > 1 requires an integration by parts 
argument similar to the one in Proposition 4.2. We consider \\xd x (e~ f — f)\\o- Using 
stochastic completeness of the heat kernel, we find 

F := xd x (e- tA f -f) = J (xd x H) ■ f(x, y, z^xT^dx dvol aM (^) 

- / {xd x )[H(t, x,y,z,x,y,z)f(x,y, z)]x~ b ~ l dx <\vo\ dM (x) =: F x - F 2 . 

Next we transform to projective coordinates and integrate by parts in S, where the 
boundary terms lie away from the diagonal and hence are infinitely vanishing for t — > 
by the asymptotic behaviour of the heat kernel. Omitting these irrelevant terms, we 
obtain 



Ft = J Ul-vS^ds-Udu + xdA [ V - m G(r],S,U,Z,x,y, 
x / (x(l - rjS),y - xrjU,z - rjZ) hr) m {l - rjSy^dS dU dZ 
= J [{-Udu +xd x )r]- m G] ■ fhrf^l-riSy^dSdUdZ 

-J G- (d s f) ■ hr]-\l - r]S)- b dS dU dZ 

-J G ■ f ■ d s r]-\h{l - r]S)- b )dS dU dZ. 

We perform similar computations for F 2 : 

F 2 = J ( K {xd x H)f{x,y,z) + H-{xd x f))x- b - 1 dxdvo\ dM {x) 

= J m - T]S)^d s - Udu + xd x ^j [r]- m G]f ■ hr] m (l - rjSy^dS dU dZ 

+ J G(7],S,U,Z,x,y,z)(xd x f(x,y,z))h(l-r ] S)- b - 1 dS dU dZ 

= J [(-Udu + xd x )r]- m G] ■ f(x,y,z)hri m {l-r ] S)- b ~ 1 dSdUdZ 

-J G- f ■d s (r)- 1 h(l-r ] S)- b ) dSdUdZ 

+ J G(ri,S,U,Z,x,y,z)(xd x f(x,y,z))h(l-r]S)- b - 1 dSdUdZ. 
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Carefully recalling the variables we have suppressed, we see that F = F\ — F 2 becomes 
F = J [(-Ud u + xd x )r ] - m G(r ) ,S,U,Z,x,y,z)]hr ) m (l-r ) S)- b - 1 
x (f(x(l — r)S),y — xrjU, z — r]Z) — f(x, y, z))dS dU dZ 
- J G^S^Z^z^ds^Kx^ - r)S),y - x V U,z - V Z)(1 - vS)-')) 
x (f(x(l - r]S),y - xrjU, z - rjZ) - f(x,y, z))dS dU dZ 
-J G-((l- r)S)d s (r]~ 1 f(x(l -r)S),y- xrjU, z - r]Z)) + xd x f(x, y, z)) ■ 

x hil-rjSy^dSdUdZ. 

All three integrals are now estimated as in the k = case. Note that in the second 
integral, the apparently singular rj^ 1 is cancelled upon first performing the differentia- 
tion; in the third integral, we use the fact that / e Cl(M,g). Higher order and other 
edge derivatives may be similarly estimated. □ 
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